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We derive nonstrange baryon-baryon scattering amplitudes in the nonrelativistic 
quark model using the "quark Born diagram" formalism. This approach describes the 
scattering as a single interaction, here the one-gluon-exchange (OGE) spin-spin term 
followed by constituent interchange, with external nonrelativistic baryon wavefunc- 
tions attached to the scattering diagrams to incorporate higher-twist wavefunction 
effects. The short-range repulsive core in the NN interaction has previously been 
attributed to this spin-spin interaction in the literature; we find that these pertur- 
bative constituent-interchange diagrams do indeed predict repulsive interactions in 
all I,S channels of the nucleon-nucleon system, and we compare our results for the 
equivalent short-range potentials to the core potentials found by other authors using 
nonperturbative methods. We also apply our perturbative techniques to the NA 
and AA systems: Some AA channels are found to have attractive core potentials 
and may accommodate "molecular" bound states near threshold. Finally we use our 



Born formalism to calculate the NN differential cross section, which we compare with 
experimental results for unpolarised proton-proton elastic scattering. We find that 
several familiar features of the experimental differential cross section are reproduced 
by our Born-order result. 
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I. INTRODUCTION 



The nucleon-nucleon interaction exhibits a strongly repulsive short- distance core and a 
longer-ranged but much weaker attraction. Although there has been evidence of the general 
features of this interaction for over fifty years the physical mechanisms proposed as the 
origin of the interaction have changed as our understanding of the strong interaction has 
progressed. In 1935 Yukawa suggested that the finite-ranged nuclear attraction was due 
to the exchange of a massive, strongly-interacting meson, by analogy with electron sharing as 
the origin of chemical forces. This hypothetical meson was identified with the pion after its 
discovery, and as the lightest hadron it certainly contributes the longest-ranged component 
of the nucleon-nucleon strong force. The repulsive short-range core of the nucleon-nucleon 
interaction was similarly attributed to the exchange of heavier mesons such as the u after 
their discovery. These meson-exchange models have been elaborated considerably since these 
original suggestions, and the most accurate phenomenological descriptions of the nucleon- 
nucleon interaction at present are meson exchange models ||, with parameters such as 
meson-nucleon couplings fitted to experiment. 

A literal attribution of the short-range repulsive core to vector meson exchange, as op- 
posed to a phenomenological parametrization, of course involves a non sequitur H||: Since 
the nucleons have radii of ~ 0.8 fm, and the range of the vector-exchange force is l/m u « 0.2 
fm, one would have to superimpose the nucleon wavefunctions to reach the appropriate in- 
ternucleon separations. The picture of distinct nucleons exchanging a physical u meson at 
such a small separation is clearly a fiction, and a realistic description of the short-range core 
interaction requires a treatment of the quark wavefunctions of the interacting nucleons and 
a Hamiltonian which involves quark and gluon degrees of freedom. 

Since the development of QCD as the theory of the strong interaction there have been 
many studies of the NN interaction in terms of quarks and gluons. Most have employed 
the nonrelativistic quark potential model, although some early work used the MIT bag 
model ||. Many studies of this and other multiquark systems were fundamentally flawed 
due to an inadequate treatment of the color degree of freedom, or due to assumptions that 
imposed confinement on the entire multiquark system; a discussion of these problems is 
given by Isgur j|. It now appears that a pairwise A • A color interaction [7|-[J together 
with a sufficiently general spatial wavefunction that allows dissociation into color singlets 
provides a sufficiently realistic description of color forces in a multiquark system. The NN 
references summarised here all assume the A • A form, and it has also become standard to 
employ a quadratic or linear confinement potential. Finally, the spin-spin hyperfine term 
is incorporated in all these references, as it apparently makes a dominant contribution to 
the NN core interaction. Several groups have included other terms from the Breit-Fermi 
interaction such as the hyperfine-tensor and spin-orbit interactions. 

Attempts to describe the NN interaction using quark potential models date from the work 
of Liberman in 1977 ||10|| , who calculated the adiabatic response of the six quark system to 
variations in the interbaryon coordinate. The resulting effective potential had a repulsive 
core with a weak intermediate-range attraction. Liberman concluded that the repulsive core 
was predominantly due to a combination of the Pauli principle and the contact hyperfine 
term. The same conclusion was reached by Neudatchin and collaborators |LT|| , also within 
the adiabatic approximation. Harvey |L2| continued this adiabatic approach in a generalized 
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calculation of the effective NN potential; he noted that SU(4) spin-isospin symmetry required 
that NN, A A, and "hidden color " (ggg)s <8> (qqq) 8 states be included in the calculation. On 
incorporating these states, he found that the repulsive core was strongly suppressed. It is 
now widely believed that the absence of the repulsive core in improved calculations is an 
artifact of the adiabatic approximation |13],[n] . 



The usual method for improving on the adiabatic approximation in the g 6 system is to 
employ the resonating-group method. This involves expanding the wavefunction of the sys- 
tem in a basis which describes system subclusters multiplied by unknown functions of the 
intercluster coordinates. One then solves the resulting coupled integro-differential equations 
using various numerical techniques. Baryon subcluster wavefunctions are usually taken to 
be simple Gaussians, and the coupled system is truncated at the NN, NN + A A, or NN 
+ AA + hidden color levels. The first applications of the resonating-group method to the 
NN system were given by Warke and Shanker ||15|| , Oka and Yazaki [16]| and Ribeiro [17 



All groups found a repulsive core, which was dominantly due to the Pauli principle and the 
color hyperfine term. The AA and hidden color channels were found to make only small 
contributions to the hard-core S-wave phase shifts. These studies found that several contri- 
butions including the Coulomb and confinement kernels approximately cancelled, leaving the 
spin-spin OGE hyperfine term as the dominant interaction. The resonating-group approach 
has been extended to include more channels, strange quarks, effective one-boson-exchange 
long-range potentials, and virtual excitations of the quark wavefunctions. In particular, 
Koike |18j has applied these techniques to the "flip-flop" model |19| (which eliminates long- 
range color van der Waals forces), supplemented by an effective meson-exchange potential, 
and Cao and Kisslinger [^0] have developed a relativised resonating-group formalism and 
applied it to the determination of equivalent potentials and low-L phase shifts in a model 
which incorporates OGE and meson exchange forces. Both these references find reasonably 
good agreement with experimental low-energy NN phase shifts. A summary of work in this 
field to 1989 has been given by Shimizu 



Maltman and Isgur [2^] have performed a detailed variational calculation of the ground 
state properties of the deuteron using a A ■ A quark-quark interaction with the full OGE 
color Breit-Fermi interaction, a quadratic confinement term, and a phenomenological one- 
pion-exchange potential. In contrast to typical resonating-group calculations, they allowed 
spatial excitations within the clusters. They found a repulsive core, and noted that the 
admixture of P-wave color octet clusters significantly increased the range and depth of the 
intermediate-range attraction. Their results for the deuteron binding energy, RMS radius, 
quadrupole moment, and magnetic moment all agreed well with experiment. 

Although the origin of the nucleon-nucleon force at the QCD level is now reasonably 
well understood, the resonating-group and variational techniques which have been employed 
in this work are rather intricate and require considerable theoretical effort, and usually 
lead to numerical rather than analytical results. As these techniques are best suited to the 
determination of ground state properties, topics such as resonance production and scattering 
cross sections at higher energies have received little attention in these quark model studies. 

In recent work we have investigated the possibility that these low-energy nonresonant 
hadronic scattering amplitudes may actually be dominated by simple perturbative processes; 
if so, it may be possible to derive useful estimates of these amplitudes using a much simpler 
approach. A complementary possibility of perturbative dominance of hadronic scattering 
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processes at high energies through constituent-interchange mechanisms has been investigated 
by theorists almost since the development of QCD p3| . Results for elastic hadron-hadron 
scattering amplitudes, in particular the asymptotic Q 2 dependence of fixed-angle scattering, 
have recently been presented by Botts and Sterman |24| . 

Concerns regarding the range of validity of the high-energy perturbative QCD stud- 
ies have been expressed by Isgur and Llewellyn-Smith |p5| , who suggest that higher-twist 
hadron wavefunction effects may actually dominate perturbative QCD contributions at ex- 
perimentally accessible energies. We have explicitly incorporated these wavefunction effects 
in our study of a constituent-interchange scattering mechanism in the nonrelativistic quark 
potential model. We calculate the hadron-hadron scattering amplitudes which follow from 
one gluon exchange followed by constituent interchange (quark line interchange is required 
at lowest order in a s to restore color singlet final states), with nonrelativistic quark model 
wavefunctions attached to the external lines. This OGE+CI mechanism may be dominant 
in processes in which qq annihilation is forbidden for the valence wavefunctions. We have 
applied this description of scattering to elastic 1=2 tttt [^6| and 1=3/2 Kir [27| reactions and 



found excellent agreement with the experimental S-wave phase shifts given standard quark 
model parameters. Related approaches to calculating meson-meson scattering amplitudes 
which iterate this quark-gluon mechanism have been discussed in the literature [p8| |. These 



Born-order techniques have also been applied to vector- vector meson systems [29|, and lead 



to interesting predictions of vector- vector molecule bound states in certain channels ||30|| . In 
the vector-vector system the hyperfine interaction apparently does not dominate the scat- 
tering amplitude, unlike the pseudoscalar-pseudoscalar and NN systems. More recently we 
applied the quark Born formalism to KN scattering ||31|| , which is also free of qq annihilation 
at the valence quark level. We found satisfactory agreement with the experimental S-wave 
KN scattering lengths, although 1=0 is not yet very well determined experimentally. The 
KN S-wave phase shifts at higher energies, however, are not well described; they require 
stronger high-momentum components in the nucleon wavefunction than are present in the 
single Gaussian forms we assumed. The higher-L KN partial waves, especially the P-waves, 
show evidence of a spin-orbit interaction which does not arise in single-channel spin-spin 
scattering, which has not yet been adequately explained in the literature. 

The next level of complexity in Hilbert space is the g 6 baryon-baryon sector. Since 
this system is free of annihilation at the valence level, and the spin-spin hyperfine term 
has already been established as the dominant interaction underlying the core repulsion, 
derivation of the NN core interaction is an important test of the quark Born formalism. Here 
we derive the nucleon-nucleon interaction predicted by the OGE spin-spin term using quark 
Born diagrams, and show that the predicted core interaction is indeed strongly repulsive in 
all four spin and isospin channels. Low energy Born-equivalent NN core potentials are also 
derived and compared to previous results. We then consider other nonstrange baryons and 
derive the NA and AA short-range interactions; some of these are found to be attractive, and 
we investigate the possibility that these channels might support dibaryon molecule bound 
states. Some of our results for attractive AA channels are consistent with the previous 



conclusions of Maltman |32]. Finally we derive the elastic NN differential cross section 
predicted by our quark Born formalism and find that some familiar experimental features 
of the high-energy elastic proton-proton differential cross section are evident in our results. 
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II. DERIVATION OF SCATTERING AMPLITUDES 



a) Hamiltonian and hadron states 

In the quark Born diagram formalism we derive the matrix element of the interaction 
Hamiltonian between quarks in incoming hadron states to leading Born order, which is then 
used to calculate scattering observables. We factor out the overall momentum conserving 
delta function and then derive the remaining matrix element, which we call tiff, 



f 



{BB'\H scat \BB')i = h fi 5(P f - Pi) 



(1) 



Since the hadron state normalizations we will introduce are identical to those used in our 
previous study of 1=2 tctt scattering [^6| we can use the relations between the scattering 
matrix element hfi and the phase shifts and cross sections given there. The details of our 
diagrammatic procedure for determining hfi are described elsewhere p6| , |27|j3T[| ; here we shall 
simply recall some basic points and then give our results. 



For baryon-baryon scattering we shall follow previous studies [p|,p!0|- pJ| , pl| ,p2| and assume 
that the dominant part of the core interaction derives from the spin-spin color hyperfine 
term, 



H 



sea t 



a,Kj 



8ira s 
3rriimj 



Si ■ Sj 



(2) 



where is the color matrix A a /2 for quark i. The baryon color wavefunctions are the usual 
color singlets, 



\baryon) = ~7f e *jk 

j,*=i,3 v6 



(3) 



Our spin-flavor states for the meson and baryon are the usual SU(6) states, but as explained 



in reference [31] we find it convenient to write these states using field theory conventions 
rather than in the usual quark model form. The quark model conventions show explicit 
exchange symmetry by assigning a fixed location in the state vector to each quark. The 
field theoretic convention greatly reduces the number of terms encountered in our scattering 
matrix elements; for example, the proton state in field theory conventions has only 2 terms 
instead of the usual 9 for quark model states, so PP— >PP elastic scattering involves only 16 
terms, far fewer than the 6561 we would encounter with the usual quark model conventions. 



As examples, the orthonormal S z 
conventions are 



3/2 |Ag/ 2 ) and 5^ = 1/2 \P\/2) states in field theory 



and 



l^i/2) 



v 3/2/ 




V2 



\u + u + d + ) 



u + u + dJ) 
~V2~ 



- \u + u^d+) 



(4) 



(5) 



The other baryon states considered in this paper can be derived from these by application 
of spin and isospin raising and lowering operators. 



6 



We shall quote general results for baryon-baryon scattering amplitudes with arbitrary 
spatial wavefunctions attached to the external lines and then specialize to single-Gaussian 
forms to derive representative closed-form results. The general spatial baryon wavefunction 
we assume is of the form 

®baryon(Pl,P2,P3] Plot) = (fibaryon (Pi , Vl , P3 ) S(P tot - fit - p 2 ~ p 3 ) (6) 

with a normalization given by 

{^baryon (-ffot ) I ^baryon (-^toi ) ) 




= 5(P tot - P;j . (?) 

The standard quark model single-Gaussian baryon wavefunction we shall use for illus- 
tration is 

, - 3 3/4 J (Pi +pi +pi -Pl-P2-P2-P3-P3-Pl)\ / a s 
<Pbaryon{pi,P2,P3) = ^3/2^3 eX P ( 3^2 J " ( 8 ) 

Oscillator parameter values of 0.25 GeV < a < 0.42 GeV have been used in the quark model 
literature on baryon spectroscopy, as we will discuss subsequently. 

b) Baryon-baryon scattering amplitudes 

By analogy with our study of KN scattering |3l| (see especially section lib) we first 
write a generic scattering diagram with initial and final baryon-baryon states. We then 
connect the initial and final quark lines in all ways consistent with flavor conservation; for 
example A ++ A ++ elastic scattering has 6!=720 quark line diagrams. These may be grouped 
into four sets in which the number of qq pairs which cross in t-channel is zero, one, two or 
three. We then generate scattering diagrams by inserting one-gluon-exchange interactions 
between all pairs of initial quarks in different initial baryons; this gives nine times as many 
scattering diagrams as we had quark line diagrams. Many of these diagrams are trivially 
zero; these include the zero-pair-interchange and three-pair-interchange diagrams, which 
vanish due to color. The nonzero scattering diagrams may be related to a small "reduced 
set" of diagrams by permutation of external lines, which leaves a diagram invariant. In 
baryon-baryon scattering this reduced set contains eight independent diagrams, which are 
shown below. 
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(10) 
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In a given baryon-baryon scattering process the matrix element hfi in (JT]) is a weighted 
sum of the eight spatial overlap integrals represented by the diagrams Di . . . D s , 



k fi = J2 w n I n (D n ) . 



(17) 



n=l 



The weight of each diagram (introduced in |!T|) is the product of a color factor, a fermion 
permutation phase called the "signature" of the diagram (which is — 1 for D\ . . . D4 and +1 
for D 5 . . .D 8 ), the overall (— ) in Hj (g), and a reaction-dependent spin-flavor factor. The 
derivation of these factors and the spatial overlap integrals they multiply is discussed in detail 
elsewhere PS|J5T| ], so here we will simply present results with minimal discussion. There is a 
minor change in the convention for diagram weights relative to our earlier reference. In our 
KN study [0 we incorporated the (— ) phase of Hi and the signature phase in the spatial 
overlap integral; here we include them in the diagram weight. This overall factor was (+1) 
for all KN diagrams, so the KN weights are unchanged by our new convention. Similarly, in 
our first paper we incorporated the (— ) phase of Hj in the spatial overlap integral. Our 
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new convention is useful because it makes all the NN spatial overlap integrals considered 
here positive, so the overall amplitude phases are clear from the weights alone. 
The color factors of the diagrams D\ . . . D s are 



Icolor([Dl ...D 8 }) 



+ 4/9, -2/9, -2/9, +1/9, +4/9, -2/9, -2/9, +1/9 



(18) 



D n is related to D n+ 4 by t <-> u crossing; for this reason the weights w n and w n+ 4 are closely 
related in many reactions. To simplify our presentation, when possible we will just give 
results for the weights of D\ . . . D4 and indicate the relative phase of the set for D§ . . . D 8 
after a bar. Thus for {I C oior(D n )} above we write 



' color 



+ 4/9, -2/9, -2/9, +1/9 (+) 



(19) 



The spin-flavor weights (incorporating the signature phases) are just matrix elements of 
the operator Si ■ Sj between two initial quarks for the given process. As an example, for 
A ++ A ++ , S — 3, S z — 3 scattering there are four \u + u + u + ) / a/6 external baryons, and only 
the SfSj terms contribute. On summing over all scattering diagrams in this channel we find 



Ispin— flavor ' I signature 



AAJ=3,S=3 



+ 9/4, +9/2, +9/2, +9 (-) 



(20) 



which is I z - and ^-independent. Combining these we find the / = 3, S = 3 AA diagram 
weights, 



{w n (AA, I = 3,S = 3)} 



1,-1,-1,+! 



(21) 



Since the spatial overlap integrals I 5 (D 5 ) . . . I 8 (D 8 ) are equal to the integrals 
ii(-Di) . . . J 4 (D 4 ) after t <-> u crossing, the relative (— ) phase in the diagram weights ( |2T| ) 
insures that the I = 3, S = 3 A A scattering amplitude hfi is spatially antisymmetric, as 
required for a totally antisymmetric fermion-fermion scattering amplitude that is symmetric 
in the remaining degrees of freedom (I and S). This antisymmetry is a nontrivial check of our 
spin-flavor combinatorics, since it is only evident after the sum over individual quark-gluon 
scattering diagrams is completed. 

The diagram weights for all I,S channels of NN, NA and AA elastic scattering are 
tabulated at the end of the paper; these and the overlap integrals constitute our central 
results. 

The spatial overlap integrals associated with the diagrams may be determined using the 
simple diagrammatic techniques presented in Appendix C of reference p6[ . In these integrals 



all momenta implicitly three-dimensional, and have an overall spin-spin coefficient n ss of 



K s 



3ml (2tt) 3 



(22) 



This is ( — 1) times the k of 0], since we have chosen to include the (— ) phase of Hi (|2|) in 



the diagram weight factor, as discussed above. The integrals are 

h — K ss J J da\d(i2 $^(ai, 0,2, A — a\ — a 2 )^Q(ai, C — ai — a 2 ) 
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db 2 db 3 ® B (-A -b 2 - b 3 , b 2 , b 3 )<5>* D (-C - b 2 - b 3 , b 2 , b 3 ) ; (23) 

h = k ss JJJJ da x da 3 db 2 dc 3 ® A (ai, A-a x - a 3 , a 3 )$^(a 1 , C - a x - c 3 , c 3 ) 

$ B (C -A + a 3 , b 2 , -C-a 3 - b 2 )<S>* D (a 3 , b 2 , -C - a 3 - b 2 ) ; (24) 
-^3 = Kss JJJJ da\da 2 db 3 ddi <& A (ai, a-2, A — a\ — a2)$^(cti, ct2, C — a\ — a 2 ) 
$ B (C - ai - a 2 , -A - C + a x + a 2 - b 3 , b 3 )<S>* D (d h -C - b 3 - d u b 3 ) ; (25) 
h = k ss / / / / da\da 3 dbidb 3 & A (ai, A — a\ — a 3 , a 3 )$^(ai, C — a\ — bx, b\) 



$ B ( 6lj _A - &! - 6 3 , & 3 )$J,(a 3 , -C-a 3 - & 3 , b 3 ) . (26) 

We evaluate these in the cm. frame, so the £ u crossed integrals I 5 . . . I 8 can be obtained 
by exchanging C and D = —C, or in terms of the cosine of the cm. scattering angle 
fx = cos(0 c . m .), 

J n+4 (/i) = I n (-/i) • (27) 

A simplification follows if all baryons have the same spatial wavefunctions, as we assume 
here; in this case I 2 = I 3 and hence Iq = I-j. 

The overlap integrals may be carried out in closed form given single-Gaussian wavefunc- 
tions (H), and each gives a result of the form 

4 = k ss i] n exp I - (A n - B n fi)P 2 ^ , (28) 

where P is the magnitude of the cm. three-momentum of each baryon, P 2 = A 2 = B 2 = 
C 2 = D 2 . The results are 

h = K ss exp| - ^-L(!-C) 2 j = K ss exp{ - -1^2(1 - /i)P 2 | = Kss exp j^j ; (29) 
!, = !. = *. (^) 3/2 exp{-^(20-12,)P 2 } 



/12\ 3/2 f 2(s-4M 2 )l f 2t I 
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In the final expression for each integral we have substituted for P 2 and /i in terms of the 
Mandelstam variables s and t using relativistic kinematics, s = 4(P 2 + M 2 ) and t = —2(1 — 

Near threshold the overlap integrals are comparable in magnitude, but at higher energies 
their behaviors differ markedly. All but I\ and its crossing-symmetric partner I5 are strongly 
suppressed in s; the diagrams D\ and D§ therefore dominate at high energies, for forward 
and backward scattering respectively. This behavior is due to the mechanism of "minimum 
spectator suppression" , as was discussed in detail in section He of reference |3l] . To summa- 
rize the arguments for this case: 1) for forward scattering, A = C, diagram Di (9) requires 
no spectator to cross into an opposite-momentum hadron, which would carry considerable 
suppression due to a small wavefunction overlap; only the hard-scattered constituents are 
required to reverse momentum. 2) D 5 (13) requires all spectators that were initially in a 
baryon with momentum A to reside finally in a baryon with momentum D. Clearly the 
suppression due to wavefunction overlaps of the spectators will be less important if the final 
baryon D has the same momentum as the initial baryon A, D = A. This corresponds to 
backscatter, C = —A, since C = —D in the cm. frame. As in Di, only the hard-scattered 
quarks are then required to recoil into a baryon with three-momentum opposite to that of 
their initial baryon. These two explanations are actually equivalent because D\ and D§ are 
related by crossing. 



III. NN CORE POTENTIALS AND PHASE SHIFTS 

The Hamiltonian matrix elements 



E W nln (32) 



n=l 



for the four I,S channels accessible in NN scattering are summarised by the diagram weights 
in Table I. Specialising to the even-L channels I,S=0,1 and 1,0, for which a repulsive core in S- 
wave is a well known feature, we see that all eight coefficients {wi . . . w^} are positive or zero 
in both cases, corresponding to a repulsive interaction. For a more quantitative evaluation, 
we can relate this matrix element to an NN potential near threshold, which is defined to 
give the same low-energy scattering amplitude near threshold in Born approximation. (See 
Appendix E of reference ]2B| for a detailed discussion.) For an hfi of the form 

87ra s 1 f ,1 

h fi = 3m 2 ( 2?r )3 ^ WnVn exp { ~ ( An ~ Bn ^> P ) ( 33 ) 

the Born-equivalent potential is 

T/ ( \ V" Wn ^ n ! ^ \ MA\ 
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The t u crossed diagrams -D5 . . . D% are not included in this sum because they will 
automatically be generated by the crossed diagram in NN— >NN potential scattering through 
V NN (r). 

The numerical potentials predicted for S-wave 1=0 and 1=1 NN systems are shown in 
Fig.l for our "reference" set of quark model parameters |3~H1 , a s = 0.6, m q = 0.33 GeV 
and a = 0.4 GeV. Actually only the two parameters a s /m 2 q = 5.51 GeV -2 and a = 0.4 
GeV are involved in V/v7v(r). These potentials are consistent with expectations for NN 
core interactions; they are repulsive and have ranges of about 1/2 fm and peak values 
comparable to +1 GeV, which is essentially infinite from a nuclear physics viewpoint. It may 
be interesting in future work to parametrize the amplitudes associated with each diagram 
(the weights in Table I) as a two-nucleon spin- and isospin- interaction of the form A I + B Si ■ 
S2 + C f\ ■ r 2 + D Si ■ S2 T\ ■ r 2 , which will allow a more direct comparison with meson-exchange 
models [03| . 

Low-energy equivalent NN core potentials have been presented as the results of some of 
the NN resonating-group and variational calculations we discussed in the introduction. In 
their Figs.l and 2 Suzuki and Hecht [Rl] show numerical results for the NN core potentials 
of Harvey fli~2 |, Faessler, Fernandez, Liibeck and Shimizu |13j and Oka and Yazaki ||T6"| . 



The Faessler et al. and Oka- Yazaki potentials are quite similar to our potentials in Fig.l, 
with values at the origin between 0.6 and 1.0 GeV and comparable ranges. Harvey finds 
potentials with somewhat longer ranges, which Suzuki and Hecht attribute to his choice of 
a larger nucleon width parameter, = 1/a = 0.8 fm; the Faessler et al. and Oka- Yazaki 
values are 0.475 fm and 0.6 fm respectively, and we use a comparable 1/a = 0.493 fm. 

The NN core potentials found by Maltman and Isgur, in Fig.l of reference [22], also have 
similar ranges but are somewhat larger in magnitude, V(0)=1.2 GeV for 1=0 and 2.3 GeV 
for 1=1. Our difference in the contact values is due in part to the choice of parameters; in 
our calculations the potentials are proportional to a s a z /m? q , which is 0.353 GeV with our 
parameters and 0.489 GeV for Maltman and Isgur. Note, however, that no other reference 
finds the large Maltman-Isgur 1=1/1=0 ratio at contact. The value chosen for a s a 3 /m 2 
by itself does not explain the differences between potentials; those reviewed by Suzuki and 
Hecht use a s a 3 /m 2 q rs 0.55 GeV, so we would naively expect our potentials to be ~ 0.6 times 
as large as theirs. Of course the values near the origin have little physical relevance due 
to their small Jacobean weight, and in any case we are comparing potentials derived using 
three different methods, and these differences may preclude a more accurate comparison of 
results. 

The choice of parameters will be discussed in more detail in the section on differential 
cross sections. Here we simply note that the smaller value of a s we use is now generally pre- 
ferred because recent spectroscopy studies using improved wavefunctions have considerably 
lowered the value required to fit hadron spectroscopy. The NN references we compare with 
predate the improved spectroscopy studies and thus used a rather large value of a s , which 
was required to give a realistic NA splitting given single-Gaussian wavefunctions. We also 
prefer to use our fixed parameter set because these values were found to give reasonable 



results for low-energy S-wave titi, Ktc and KN scattering in our previous studies |26| , ^7| , |31 

Note in Fig.l that the intermediate-range attractions which are responsible for the 
deuteron in 1=0 and its almost-bound 1=1 partner are absent from our quark Born po- 
tentials. This is as we anticipated, given that these attractions arise mainly from a spatial 
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distortion of interacting- nucleoli wavefunctions |4,[22|; in our leading-order Born calculation 



we assume fixed nucleon spatial wavefunctions. The attraction presumably arises at higher 
order in the Born series, and may be accessible through leading-order Born calculations of 
off-diagonal matrix elements. 



Oka and Yazaki (Fig.2 of ||16| ) and Koike (Fig. 3 of fTS| ) also show the S-wave phase shifts 
which result from their NN core interactions. In Fig.2 we show the S-wave phase shifts we 
find on numerically integrating the Schrodinger equation with the potentials of Fig.l. Our 
phase shifts are very similar to the results of these earlier resonating-group studies. Although 
we would like to compare our phase shifts to experiment directly, the experimental phase 



shifts |3^J are unfortunately complicated by the presence of the deuteron and its 1=1 partner 
near threshold. These states will have to be incorporated in our calculation before we can 
make a useful comparison between our theoretical core phase shifts and experiment. 

The Born-order approximate phase shifts (proportional to hfi) can be determined ana- 
lytically using Eq.(6) of |2TJ and dividing by 2 for identical particles, 

S® = -^PE P £ h fi (ji)P t (ji). (35) 

The momentum, energy and fi = cos(# c . m .) are for one nucleon in the cm. frame. From 
our general result for hfi ( |P7|) and the Gaussian- wavefunction integrals fl28|) , we find an £ih 
Born-order partial-wave phase shift of 

g 

*W = -^PE P £ w nVn e- A " p \ e (B n P 2 ) . (36) 

6m q n=l 

The spin-dependence of this result is implicit in the weights {w n }. Note that these phase 
shifts are functions of i and spins only, so there is no spin-orbit force in our effective NN 
interaction. This is as expected given that our only interaction at the quark level is the 
spin-spin hyperfine term. A more realistic model will require a generalization to include the 
OGE spin-orbit term and perhaps coupled channel effects, as we discussed in our study of 
KN scattering. 

The analytic Born-order result (|36"D for the phase shifts is unfortunately of little utility 
for NN S-waves given realistic quark-model parameters; the equivalent potential V^^^r) is 
nonperturbatively large in this case, and must be iterated coherently to determine phase 
shifts, as we have done in Fig.2 using the NN Schrodinger equation. In contrast to this 
result, we previously found nonperturbative effects in the S-wave phase shifts of tttc, KK 
and KN systems to be much less important, due to somewhat shorter-ranged forces and 
the smaller reduced mass. The NN Born-approximation phase shifts (|36|) are presumably 
more useful for higher partial waves and higher energies, since multiple scattering effects are 
expected to be largest in S-wave near threshold. 



IV. OTHER NONSTRANGE BB' CHANNELS: NA AND AA 

The core potentials predicted by the quark Born formalism for other nonstrange baryon- 
baryon channels should allow tests of the assumed hyperfine dominance in systems other 
than the familiar S-wave NN cases. The short-range interactions in the NA and AA channels 
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may be observable experimentally as final state interactions or, if the interaction is suffi- 
ciently strong to support bound states, as dibaryon molecules not far below threshold. The 
possibility of nonstrange resonances in the g 6 sector has been considered by many authors, 



the earliest reference apparently being a group-theoretic study by Dyson and Xuong p6 
These q 6 systems have also been studied using the bag model (which is unfortunately known 
to give unphysical predictions of a host of multiquark resonances) , one-boson-exchange mod- 
els and the nonrelativistic quark model; references before 1985 are summarised by Maltman 
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Our results for NA and AA are summarised by the diagram weights in Tables II and 
III. After completing their derivation we found that some of these matrix elements had 
previously been tabulated by Suzuki and Hecht p4| ; our NN and AA weights wi, u>2, W3, W4 



are equivalent to the coefficients C^j., — C^j., — Cjj., in their Table II, which provides an 
independent check of our results in these cases. As these weights multiply comparable spatial 
overlap integrals which give positive contributions to the low-energy equivalent baryon- 
baryon potential fl3"4|), negative weights imply attractive potential contributions. 

Referring to Table II, we see that the NA system has a strongly repulsive core in the 
channels I,S=2,1 and 1,2 and weak core interactions in 2,2 and 1,1. Unlike the AA system 
(to be discussed subsequently) our NA core interactions do not lead to bound states in any 
channel. Since the lightest reported dibaryons have masses very close to the NA threshold 
35| , |37|j38| , the experiments may be seeing threshold effects due to the opening of the NA 



channel, or perhaps weakly-bound NA molecules. Our calculation does not support the 
existence of such bound states, although the NN system is similarly predicted to have a 
purely repulsive core, but the I,S=0,1 deuteron is nonetheless bound by an intermediate- 
range attraction which is absent from our leading-order Born calculation. Similar weakly- 
bound states may exist in NA and AA as well, despite repulsive cores. 

Next we consider the AA system. For AA some general rules follow from the assumption 
of a single A • A interaction; since the initial three-quark clusters are transformed into color 
octets by the interaction, line diagrams with zero or three pairs of quarks exchanged are 
forbidden. Thus the amplitude for A ++ A~ elastic scattering must be zero. This implies 
relations between AA amplitudes with different isospins, 

h fi (AA; I = 0,S) = -h fi (AA; I = 2,S) (37) 

and 

h fi (AA; 1 = hS) = -I MAA; I = 3, S) . (38) 

9 

Specializing to the (+)-symmetry (even-L) cases in Table III, which include the S-wave 
channels that are a priori the most likely to support bound states, it is evident that two 
AA channels have strongly attractive core potentials, I,S=1,0 and 0,1. Of these 0,1 has 
the strongest attraction. We search for bound states by solving the Schrodinger equation 
in the A A system using the low-energy potential (|34]). With our reference parameter set 
a s /m 2 q = 5.51 GeV -2 and a = 0.4 GeV the attractive core is too weak to induce binding. 
Note, however, that previous studies of baryons using single-Gaussian wavefunctions have 
generally assumed a much stronger hyperfine term, for reasons we will discuss subsequently. 
If we use a typical parameter set from these references, a s /m 2 q = 14.9 GeV -2 and a = 0.32 
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GeV (Maltman and Isgur p2[), we find a single S-wave A A bound state in the 0,1 channel, 
with Eg = 40 MeV. Of course this channel has a fall- apart coupling to NN, so a coupled- 
channel treatment including the NN system may be required to search for resonant effects. 
None of the other AA channels have sufficiently strong attractive cores to form bound states 
in our formalism with the Maltman-Isgur parameters. 

Our result for the I,S=0,1 AA channel is remarkably similar to the conclusion of Maltman 
52 1, who found that the 0,1 channel has the strongest diagonal attraction in the AA Hilbert 



space, and that these diagonal forces led to a AA bound state with E B = 30 MeV. Maltman 
concluded, however, that off-diagonal effects due to the excitation of hidden-color states 
eliminated this bound state and led to binding in I,S=3,0 (E B = 30 MeV) and 0,3 (E B = 260 
MeV) instead. The 3,0 and 0,3 channels had previously been suggested as possibilities for 
A A bound states |T6|,|39|,|4*0[] . In contrast we find strong repulsion in the 3,0 channel and a 



weak core in 0,3. 

In view of the parameter- and approximation-dependence of predictions of AA bound 
states and the theoretical uncertainties in treating hidden-color basis states, the possibility 
of nonstrange dibaryon molecules should be regarded as an open question for experimental 
investigation. The channels which appear of greatest interest at present are the attractive- 
core systems I,S=0,1 and 1,0 and the 3,0 and 0,3 channels, which previous studies suggested 
as possibilities for binding. 



V. NN DIFFERENTIAL CROSS SECTIONS 



We can use Eq.28 of reference 



da 4ir 5 s 
~dt = (s - AM 2 ) 



l fi\ 



(39) 



to determine the nucleon-nucleon differential cross section in leading Born approximation, 
given the NN hfi matrix element ([H])- For the experimentally well-determined case of 
unpolarised PP elastic scattering, we have a weighted sum of the S=0 and S=l differential 
cross sections, 



da 



PP, unpolarised 



1 da 
Tdl 



^3 da 
i=i,s=o 4 dt 



7=1,5=1 



(40) 



To obtain the S=0 and S=l cross sections one simple substitutes the appropriate 1=1 dia- 



gram weights {w n } from Table I and the integrals {I n } from (|29| - |3T1) . 

Before we discuss our prediction for this differential cross section we briefly recall the 
experimental unpolarised PP result. This is shown for a range of Pi a b in Fig. 3, adapted from 
Ryan et al. [flT], Ankenbrandt et al. |42j, Clyde et al. |43|, Allaby et al. and from the 
ISR data of Nagy et al. \ 45] and Breakstone et al. |46| . The data in the figure were obtained 
from the Durham-Rutherford HEP data archive. Near threshold the angular distribution 
is approximately isotropic, but as Pi a b increases the scattering at large angles falls rapidly, 
and at high energies the differential cross section is dominated by an asymptotic "diffractive 
peak" , 
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lim 

s^oo, \t\/s«l 



da 
~dt 



cxpt. 



ae 



(41) 



For the purely hadronic part (as distinct from the divergent forward Coulomb peak) one 
finds 



and 



i expt - « 70 mb GeV -2 



(42) 



& cxpt. ^ n _ 12 Ge y- 



(43) 



for the asymptotic form [IB . 

On evaluating ( [i~0|) for NN scattering using (|17|), (|39D and Table I, we find that several 
features of the experimental differential cross section are successfully reproduced by our 
Born-order calculation. The theoretical Born-order cross section (f4"0"|) which follows from 
our "reference" parameter set a s /m? q = 5.51 GeV -2 and a = 0.4 GeV is shown in Fig. 4 for 
-Piab = 1.05,1.75,3. and 10. GeV, selected for comparison with Fig.3. Although Pj a b = 10. 
GeV superficially appears to be very relativistic, in the cm. frame it actually corresponds 
to P c .m. = 2.07 GeV, which for nucleon-nucleon scattering is only quasirelativistic. 

First note that the smooth evolution from an isotropic angular distribution to an asymp- 
totic forward-peaked one with increasing s is a simple consequence of the suppression in 
s of all diagrams at small \t\ except D\. In our calculation the contributions of the other 
diagrams fall exponentially with s. The experimental large-angle scattering does not fall 
this rapidly, and the discrepancy is probably due to our use of single- Gaussian forms; the 
actual proton wavefunction has short- distance quark-quark correlations, which presumably 
lead to power-law contributions at large s and \t\. 

Second, the observed approximate asymptotic form ([41]) is actually predicted by our 
single-Gaussian Born calculation. The overall normalization a and slope parameter b for 
this process are predicted to be 




6364 tt« 2 
19683 m\ 



(44) 



and 



3a 2 



(45) 



The theoretical result ( f44|) for the magnitude of the forward peak is actually independent of 
the spatial wavefunction, since the defining integral (|23"D is just the product of two normal- 
ization integrals in the limit A = C. 

With the reference parameter set we predict a somewhat smaller, broader peak than is 
observed experimentally, with 



12. mb GeV" 



(46) 



and 
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b = 4.2 GeV" 2 



(47) 



Both a and 6, however, are sensitive to the choice of quark model parameters, and vary by 
factors of about 10 and 3 respectively when a s /m 2 q and a are varied through a plausible 
range, which we shall discuss below. If we use typical ISR experimental intercept and slope 
values of a = 70. mb GeV -2 and b — 11. GeV -2 ^fj as input to fix our two parameters, the 
fitted values are 

^ = 13.3 GeV" 2 (48) 
m 2 q 

corresponding to m q = 0.21 GeV if we leave a s = 0.6, and 

a = 0.246 GeV . (49) 

Although these fitted parameters give the observed intercept a and slope b, the higher- 1 i | 
wings of the resulting distribution fall much too rapidly with s. This is prob ably an artifact 
of our use of soft single- Gaussian wavefunctions, and a calculation of the differential cross 
sections which follow from more realistic wavefunctions will be a very interesting exercise. 

The fitted strength of the hyperfine interaction fl48|) is now believed to be rather large, al- 
though it is similar to the values used in many previous quark-model studies of baryons. Ex- 
amples of previous values in chronological order are a s /m q = 15.5 GeV -2 (Oka and Yazaki, 
1980 [PJ]); 37.8 GeV -2 (Harvey, 1981 [13]; this value now appears exceptionally large); 7.7 



GeV" 2 (Faessler et al, 1982 0); 14.9 GeV" 2 (Maltman and Isgur, 1984 [[22J); and 14.4 
GeV" 2 (Koike, 1986 ]IS|). Large values were required to fit the N-A mass splitting given 
single- Gaussian wavefunctions; since this is proportional to (a s /m 2 )|?/>(0)| 2 , an underesti- 
mated wavefunction at contact must be compensated for by a large (a s /m q ). If one instead 
uses the actual Coulomb plus linear wavefunctions from the nonrelativistic Schrodinger equa- 
tion with m q ~ 0.3 GeV, the larger value of ^(O)! leads to a much smaller a s /m q ~ 5 GeV -2 . 
Since our scattering calculation uses Gaussian wavefunctions, one could argue which param- 
eter value is most appropriate; an improved calculation with more realistic wavefunctions 
will probably be required to eliminate these parameter uncertainties. 

The baryon width parameter a has also been assigned a rather large range of values 
in previous work. Representative values in chronological order are a = 0.41 GeV (Copley, 
Karl and Obryk, baryon photocouplings, 1969 |47|]); 0.32 GeV (Isgur and Karl, baryon 



spectroscopy, 1979 f48[); 0.41 GeV (Koniuk and Isgur, baryon photocouplings, 1980 



0.33 GeV (Oka and Yazaki, NN interactions, 1980 |§); 0.25 GeV (Harvey, NN interactions, 
1981 H); 0.42 GeV (Faessler et al, NN interactions, 1982 pf); 0.25 GeV (Hayne and 
Isgur relativised quark model, 1982 ||50||); 0.32 GeV (Maltman and Isgur, NN interactions, 



1984 |2f); 0.34 GeV (Koike, NN interactions, 1986 pf); 0.3 and 0.42 GeV, with the smaller 



value preferred (Li and Close, baryon electroproduction, 1990 [5T|j). 

Most potential models assume a value of m q near 0.3 GeV, although the relativised 
models of Hayne-Isgur fl50l , Godfrey-Isgur |52| (mesons) and Capstick-Isgur (baryons) [53| 
use a lower value of 0.22 GeV. Although these relativised models also use a small value of 
&s{Q 2 = 0) = 0.6 for the infrared limit of an effective running a s (Q 2 ), the hyperfine strength 
a s (0)/m q = 12.2 GeV -2 is again large because of the smaller m q . Its effects are reduced, 
however, by the use of a "smeared" contact interaction. 
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The m q and a used by Hayne and Isgur are essentially identical to the values we need to fit 
the slope and intercept of the experimental PP diffractive peak, although this is presumably 
fortuitous agreement. 

Although we can fit the magnitude and t-dependence of the small- \t\ differential cross 
section at high energies reasonably well with our quark Born results, we emphasize that this 
is at best an incomplete description of diffractive scattering, because the Born amplitude is 
purely real whereas the experimental small- \t\ amplitude is known to be close to imaginary 



J54J. This may imply that the first Born approximation is inadequate for small- \t\ diffrac- 
tive scattering, and that the coupling to inelastic channels is an essential component of a 
description of the diffractive amplitude, even for elastic processes |[55|| . It may be necessary 
to iterate the effect of diagram D\ to generate the observed phase |56f , perhaps including a 
sum over virtual inelastic channels. Our conclusion that suppression due to the spectator 
lines is the dominant origin of the observed diffractive t-dependence |3l| would presumably 
be unchanged by iteration of the qq hard scattering process in diagram D\. 

In view of the complexity of the baryon-baryon scattering problem, which involves a 
sum of thousands of diagrams and the evaluation of 36-dimensional overlap integrals, and 
the questionable accuracy of our nonrelativistic single- Gaussian wavef unctions, we find our 
approximate agreement with experiment encouraging. The most important discrepancies 
are in the phase of the scattering amplitude (which may require a higher-order Born study) 
and in the higher- 1 1 1 "wings" of the distribution, which may be more accurately described 
by more realistic Coulomb plus linear baryon wavefunctions. A determination of the proton- 
proton differential cross section given more realistic nucleon wavefunctions would be a very 
interesting future application of this formalism. 



VI. SUMMARY AND CONCLUSIONS 

We have applied the quark Born diagram formalism to nonstrange baryon-baryon elastic 
scattering. In this approach the hadron-hadron scattering amplitude is taken to be the 
sum of all single quark-pair interactions followed by all allowed quark interchanges, with 
nonrelativistic quark model wavefunctions attached to the external lines. This may be 
a useful description of reactions which are free of qq annihilation. The model has few 
parameters (here only two, the baryon oscillator parameter a and the hyperfine strength 
a s /rriq, since we incorporate only the OGE spin-spin hyperfine term in this study), and 
with Gaussian wavefunctions and a contact interaction the scattering amplitudes can be 
derived analytically. The model was previously applied to 1=2 nn and 1=3/2 Kn scattering 
with good results, and also gives reasonable results for low-energy S-wave KN scattering, 
although there are discrepancies at higher energies and in higher partial waves. 

NN scattering is an important test of this approach because it is also annihilation-free 
(at the valence quark level), and the baryon wavefunction and the dominance of the spin- 
spin OGE hyperfine interaction in NN are already reasonably well established. We find that 
the quark Born diagrams predict repulsive core interactions in both S-wave NN channels, 
and the equivalent low-energy potentials we extract from the scattering amplitudes are very 
similar to the results of previous resonating-group and variational calculations. We also give 
results for the NA and AA core interactions induced by the OGE spin-spin term, and find 
that certain AA channels have attractive cores and may possess bound states. Finally we 
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determine the NN differential cross section predicted by our Born amplitude, and compare 
the results with the experimental unpolarised PP differential cross section over a wide range 
of energies. We find that several well known features of experimental PP scattering are 
evident in our Born results, including the development of a high-energy forward peak with 
an approximately correct width and magnitude. 
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TABLES 



Table I. Diagram weights for NN elastic scattering. 
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Table II. Diagram weights for NA elastic scattering. 
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Table III. Diagram weights for AA elastic scattering. 



uj\ U2 ^3 <^4 rel. phase W5 . . . ujs 



1=3; S=3 


1 


-1 


-1 


1 


(") 


S=2 


5 
3 


1 

3 


1 

3 


-1 


(+) 


S=l 


19 
9 


1 

9 


1 
9 


1 
9 


(-) 


s=o 


7 
3 


1 

3 


1 

3 


1 


(+) 


1=2; S=3 


1 

3 


1 

3 


1 

3 


1 
3 


(+) 


S=2 


5 
9 


1 

9 


1 

9 


1 
3 


(-) 


S=l 


19 
27 


1 

27 


1 

27 


1 

27 


(+) 


s=o 


7 
9 


1 

9 


1 
9 


1 

3 


(_) 


1=1; S=3 


1 

9 


1 

9 


1 

9 


1 

9 


(") 


S=2 


5 

27 


1 

27 


1 

27 


1 
9 


(+) 


S=l 


19 
81 


1 

81 


1 

81 


1 

81 


(-) 


s=o 


7 

27 


1 

27 


1 

27 


1 
9 


(+) 


1=0; S=3 


1 

3 


1 

3 


1 
3 


1 
3 


(+) 


S=2 


5 
9 


1 

9 


1 
9 


1 

3 




S=l 


19 
27 


1 

27 


1 

27 


1 

27 


(+) 


s=o 


7 
9 


1 

9 


1 

9 


1 
3 
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